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^ ! Abstract 

m 

(N 

We study quantum electrodynamics coupled to the matter field on singular back- 
ground, which we call defect. For defect on the infinite plane we calculated the 
fermion propagator and mean electromagnetic field. We show that at large distances 
from the defect plane, the electromagnetic field is constant what is in agreement 
with the classical results. The quantum corrections determining the field near the 
O !* plane are calculated in the leading order of perturbation theory. 



1 Introduction 

It is a long-standing problem to formulate a quantum gauge theory with boundaries [1,2, 
3]. There are different approaches for its solution. One can impose boundary condition 
on the gauge and ghost fields [4]. Other possibility in quantum electrodynamics (QED) 
is to put "physical" boundary conditions on the electromagnetic field strength and use 
either explicitly gauge-invariant methods (see for example [5]), or consider these boundary 
condition as constrains quantizing the potential A^. From the formal point of view, 
shortcoming of the models of such a kind is that the gauge invariance, renormalizability 
or locality appear to be broken. There is another reason to doubt that these models 
describe correctly the physics of interaction of quantum field with the boundary, since 
the same sharp condition is imposed on all the modes of the fields. In reality the modes 
of the field with high enough frequencies are not constrained by the matter and are not 
influenced by its presence. We restrict ourself with these short comments about the 
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method of quantization with boundary constrains on the fields because it is not used in 
this paper 1 . 

The other approach which we follow here is to model boundary effects with a static 
background (defect) coupled to the quantum fields. The simplest background is the 
singular one. The first attempt to build renormalizable local QFT with boundaries in 
this way has been made by Symanzik in [7]. It was shown that some simple boundary 
conditions (namely Dirichlet and Neumann ones) are implemented by adding a singular 
interaction to the action density. The technique was generalized to manifolds with curved 
boundaries in subsequent work [8]. Recently the models of the scalar field theory in 2- and 
3-dimensions were considered in [9, 10, 11] (in the later paper structure of divergencies and 
corresponding renormalization procedure are also discussed). The method was successfully 
applied for the investigation of the critical behavior of magnets and alloys with free 
surfaces [12, 13]. In context of the Casimir problem the models of interaction of the 
free quantum field with singular background (S function-like potential) were considered 
by many authors [14, 15, 16, 17, 18, 19]. Many important results are known also for 
non-relativistic quantum mechanic systems with (^-potential [20, 21]. 

The problem of renormalization of QED with singular background ware considered in 
[22, 23]. The radiative QED correction to the Casimir energy was calculated in [24]. A 
review on Casimir effect was published recently [25], where most results and controversies 
in the field are discussed. 

Theory of Casimir effect based on assumption that it is a macroscopic phenomena 
generating by vacuum fluctuations of QED fields is in a good agreement with experi- 
mental data [26, 27, 28]. It is naturally to expect that the Casimir effect is not unique 
macro-manifestation of quantum fields' fluctuations, and in many situations the descrip- 
tion of macro-system behavior obtained by classical electrodynamics (ED) needs essential 
quantum corrections. The calculation of them is a theoretical task being of practical im- 
portance for development of nano-technologies and design of microdevices [29] . A possible 
method of its solution is proposed in this paper. 

A typical ED statement of problem is to find the electric and magnetic fields for given 
boundary conditions, charge and current distributions [30, 31]. In this paper we consider 
the QED version of simplest problem of such a kind. We study the gauge invariant, 
local, renormalizable model of simple defect on the plane in QED. It is suggested for 
calculation of quantum corrections for the fields of a plane in classical electrodynamics. It 
is essential that for renormalizability of the model the direct interactions of the boundary 
both with the photon and with the Dirac's fields are necessary. We calculate the leading 
order approximations for mean strengths of electric and magnetic fields expressed by 
usual relations in terms of components of electromagnetic field tensor. The results look 
like classical ones on the large distance from the plane. For small distances r — > the 
strength of fields appear to be singular as Ci/r 2 + C 2 /r, where Ci, C 2 are constants. It 
is essentially non-classical effect generated by interaction of the Dirac fields with defect 

1 Many results of application of this method for investigation of Casimir effect are presented in [6, 25] 
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on the plane. In the framework of the free scalar field theory with singular background a 
similar phenomena was studied in [32]. 

2 Statement of the problem 

We consider the QED with homogeneous defect on the infinite plane x% = invariant in 
respect to coordinate reflections. It is specified by the action functional of the form 

S de 0, ^ A; A, q, /) = S X S, + S t (A), (1) 

S Xq {ip,ip)= J ip{x,0){\ + q)ifj{x,0)dx, Si(A) = JlA(x,0)dx + J l'd$A(x, 0)dx, 

where ip, are the Dirac's spinor fields, A is the electromagnetic vector-potential, q, /, /' 
are fixed 4- vectors, q = q^ (7 M are the Dirac's gamma-matrices), and we used the short 
hand notations for the 4-vector x: x = (xo,xi,X2,xz) = (x, X3). The notations of this 
kind will also be used later. For gauge invariance of the defect action it is necessary to 
set h = 1' 3 = 2 . With this restriction the functional (1) is a most general form of local 
defect action without parameters of negative dimension and being invariant in respect to 
coordinate reflections and gauge transformations. Therefore in virtue of usual criteria of 
renormalization theory the model defined by addition to the QED action the defect action 
(1) is renormalizable [7]. It remains to be renormalizable for /' — too, because as we 
see below, there are no divergencies which need for their cancellation the /'-term in the 
Si term. 

The physical meaning of the vector I = (/, 0) is very simple. It defines the classical 4- 
current on the plane defect. By neglecting in our model the interaction of the photon and 
Dirac fields, the mean electromagnetic field coincides with solution of Maxwell equations 
with the current (supported on the plane) defined by Si. Vector q and scalar A describe 
the interaction of current and density of Dirac field with material defect. Interaction of 
vacuum fluctuations of the Dirac field with the background generates quantum corrections 
to usual classical effects. 

We calculate the leading approximation for electromagnetic field generated by the 
defect. Only the /-term in Si appears to be necessary for cancellation of ultra-violet 
divergences in our results. We set /' = and do not consider a trivial contribution to the 

— * — * 

first order effects from /'-term in Si. We choose the vector / proportional to q: I = q£, 
and show that this "minimal" form of Si with only one extra to Sf parameter £ provides 
the cancellation of divergencies by renormalization of £. 

The full action of our model has the form 

S(i>,^,A) = S QE D (jj, ij, A) + S de S, rl>, A; A, q, £), (2) 

2 We consider the gauge transformations — > + d^<j) n °t changing asymptotic of the field A^x) 
for large x. It follows from this assumption that lim Xi ^± 00 4>{x) — <fio, where 0o is a constant being the 
same for all i = 0, 1, 2, 3. Therefore Si(A) is gauge invariant. 
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where Sqed(^,iP,A) is the usual QED action: 

S QED $,4>,A) = J $(x)(id - eA(x) - m)^{x)dx f F lxu (x)F^(x)dx. 

Here, 

Fp,v = d^A u — d u A^, A = A^^, 

A^ is the potential of electromagnetic fields and gamma matrices fulfill the commutation 
relations 

with metric tensor = diag{l, — 1, — 1, — 1}. The considered model is invariant in 
respect to translations and Lorenz transformations leaving unchanged the coordinate x 3 , 
if the vectors q transformed correspondingly. 

In virtue of Lorenz-invariance of the model in the sense mentioned above, we can 
classify the defect properties with the value of the invariant q 2 = q 2 — q\ — q\. So, we have 
3 cases: q 2 = k 2 > 0, q 2 = —k 2 < and q 2 = 0, and we can choose (and restrict ourselves 
with) three coordinate systems where q — (k, 0, 0), q — (0, k, 0) q — (k, k, 0) accordingly. 

This means that interaction of QED fields with the defect is described by four param- 
eters: k, t = q 3 , A and £. 

In this paper we calculate the mean tensor F^ u of electromagnetic field: 

F, V = CJ F^e^^DAD^, (3) 

where the constant C is defined as follows 

C- 1 = J e iS( ^'^ A) DAD^Dij. 

For calculations, it is more convenient to use the Euclidean version of the model. It is 
obtained by the usual Wick rotation 

t -> -it, A -> iA , 70 -> ijo, q -> ^go- 
Then the Euclidean action can be written as 

qE qE i qE 

D — ^QED "+" ^defi 

where 

Sqed = $(id ~eA + m)^ + -F^F^, S% f $, A; A, q, £) = *(A + + £(<fl), 

= ^(z, 0), ^(x) = 1p(x, 0). 

The matrix a corresponding to 4- vector a is defined in Euclidean theory as a = X^ =0 a ^l^- 
For convenience of calculations we shall use the projection operator onto the plane x 3 = 

n(x,y) = 5(x-y)5(y3) (4) 
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and present the function *&(x), ty(x) in the form 

= J Hy)n T (y, x)dy = #} T , = / ^(x, y)iP(y)dy ee ^. 

The tensor F^ u is gauge invariant, therefore is independent of the choose of the gauge. 
We provide calculations in the Feynman's gauge 3 , using the formula 

6„ 



D^(x,y) 



J [IV 



4tt 2 (x — y)' 



for the photon propagator in configuration space. Integrating by parts in functional 
integral (3) we obtain 

•VW -^J d y (^7)4 ' w 

where 

My) = ejM - l»5{y 3 ), Uy) = C [ e- s ^^( y )^( y )Di;Dip. (6) 



In virtue of invariance of the action (2) in respect to translation of coordinates xo, x±,X2 
and reflection of x 3 , j^y) is an even function of coordinate y 3 only: J(y) = J(y 3 ) = 
J(— y 3 ), and after integration over y in (5) we obtain the following result 



1 [ x -> 

■7>0*0 = - [SusJ^ys) - S,aJ v (y 3 )]dy 3 = (7) 

Z J —x->, 



l r xs 

-x 3 

-|*a| 



= Sign(x 3 ) |e jf [<Wm(2/s) - 5^j v {y 3 )]dy 3 - £{8^% - <W*,) j , 

where Sign(xs) is the signum-function: Sign^) = X3/IX3I. 

The vector in (6) is the current generated by vacuum fluctuations of Dirac fields. In 
virtue of Farri's theorem, it vanishes in absence of defect. By usual methods of renormal- 
ization theory and its modification for the quantum field theory with singular background 
[7] one can prove that in the framework of renormalized perturbation theory jn(y 3 ) is 
presented by a sum of diagrams with all the necessary subtractions. Therefore it is finite 
(for the leading approximation it will be clear from the evident formula below), and in 
calculation of (7) there is the only problem - one with non-integrable singularity of jn{y 3 ) 
at y 3 = 0. Therefore the integral (7) needs a regularization. Since for x 3 7^ the deriva- 
tive of J-fj, u (x 3 ) in respect to x 3 is finite, we obtain the finite value of J-^x^) subtracting 
from integral in the left hand of (7) a constant dependent on the chosen regularization. 
This subtraction can be generated by term S\ with appropriate choice of the parameter 
£. Therefore the /'-term in Si is not necessary for renormalizability of considered model. 

3 This gauge can be fixed by the usual Faddeev-Popov trick 
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The leading approximation for can be presented as 

j^=Tr(S(y,y)^). (8) 
Here S(x,y) is the free fermion propagator in the theory with defect S\ q {^i,ili) 

S a p{x, y)=C f J e- s f^%(y)ip a (xWD^. (9) 
We used in (9) the notation Sf(ip, ip, q) for the fermion part of the action S E 

S f @, i/>, q) = S E ($, ij, A; A, q, £) U=o = W + *(A + K = id + m, 

and 



CJ 1 = J e- s ^ q) D^)Di). 



We restrict ourselves with calculation of the main approximation for J 7 ^, and in virtue 
of (8), the most nontrivial part of this problem is to find an evident form for the fermion 
propagator S defined by (9). 



3 Calculation of S 

The right hand side of (9) can be obtained by differentiation of generating functional 

G(fj,r)) = C f J exp{-S f $,frq)+ik + fM}v$Vil> 

where fj, rj are the fermion sources. Obviously, 

5 5 

Sofi&y) = Sfj a (x) (fy^y) ^' 7 ^'^ ' G ^ ,f ^ = exp {77^77} . 
To calculate G(fj, 77) we present the contribution of the defect in the following form 

exp{-*(A + g)*} = cyexp{CC + C(A + g)* + fC}^C, ^ = /exp{cc}^C, 
and write G(fj, rf) as 

G(r h n) = cc, J ex P {-m + CC + *, + * + *C + C(A + «*} WWW 

= cC f J exp {-$Kijj + CC + # + (hj>} V$Vi/jV(V(, 

where 

= f] + c(A + q)n, 9 = 7] + n T (. 



Now, integrating over the fields ip, ip and then over (, ( we obtain 

G(r),fj) = exp {77^77} , 

where 

S = D + S def , S def = -Dn T Q-\\ + q){lD. (10) 
We used here the following notations 

d = k~\ Q = 1 + (X + q)nt)n T . 

We see that the propagator S in our model is the sum of propagator D of usual QED and 
Sdef generated by the defect. The propagator D does not make contribution into (8), and 
we can present as 

jM(y)=Tr(S def (y,y)^). (11) 

Now, we calculate S de f in an evident form, which will be used to obtain the final result 
for T^y in considered approximation. In virtue of translation invariance for coordinates x, 
it is natural in our model to use for calculations the 3-dimensional Fourier transformation 
defined as follows 

F(x,x 3 ) = J dpe^F(p,x 3 ). 

For the propagator D we have 

D( X ) S b(i, *») = / '" + P + f/ iSnfa) exp { -E\ X3 \ + m] dp, 

where 



p = p7, E = yp 2 + m 2 . 

Hence, 

D(p,x 3 ) = ^ + P + ^73Sign(x 3 ) e _ g|g3 | 
For calculation of S de f(x,y) we note that in virtue of (4) 

nDn T (p)=D(p,0) = ^. 

So, we can write 

S def (x, y) = / ^ e ^-^^ e/ (p; x 3 , 1/3), (12) 

where _^ 

5 de/ (p; x s , y 3 ) = ~D(p; x s ) (l + (X + <7) ! ^) (A + -y s ) 



2E 2 [E (4 + q 2 + A 2 ) +4 (Am 



e -A'(|X3| + |l/3|)_ 



Here, we have denoted 

Pi = (P, iSign(x 3 )E), p 2 = (p, -iSiga(y 3 )E), 
M = 2E\ + (A 2 + q 2 )m, r = (2Eq - (A 2 + q 2 )p, 2Eq 3 ) . 

The Dirac field propagator S for the defect with A = qi = q 2 = q 3 = was calculated 
in [15, 16], and our result (13) generalize one obtained there, and coincides with it for 
particular values of defect parameters. 

Now, we use (13) for calculation of electromagnetic field generated by the defect. 



4 Calculations of J 7 , 



[lis 



Setting the right hand side of (12) into (11) and using (13), one obtains 

/°>(s) = (j {0) (x)ji°\x)) = j^-J dpTrU^ 3 )S def (p,x 3 ,x 3 )} 



(2tt)* 



-2E\x 3 \ 



{p(pq-m\)-E 2 q, o). 



E(E (4 + q 2 + A 2 ) + 4 (Am - pq)) 
After integration over angular variables j^\x) is written in the form 

j <f,) (x) = etf,0)<Kx 3 ), 

E + ap + f3m 



0(X 3 ) = . / Pi 

2ir \q\ot Jo 



-2E\x 3 \ 



2pa - (E(l - a 2 )+mf3)) In 



E — ap + f3m 



where 



a 



A\q\ 



A + q 2 + \ 2 ' 



P 



4A 



4 + g 2 + A 2 ' 



dp, (14) 



(15) 



For x 3 7^ the current j^(x 3 ) is finite since the integral (14) converges. However we 
can not set j(y 3 ) = j^(y3) in (7) directly, since j^(y 3 ) has a non-integrable singularity 
at y 3 = 0. The regularized form of the leading approximation of can be obtained 
by introducing a cut-off parameter A in (14) and substituting then of regularized current 
j(°)(A) into (7). In doing so one obtains the following result 



•F \xv\X) 



e{q^ u3 - q u 5 3fl )Sign(x 3 ) 



Hx 3 , A)- | 



where 



F(x 3 ,A) = 



1 



47T 2 |g|o; 2 JO 



A (l 

dp— 



-2E\x 3 \^ p 



E 



2pa - (E(l - a 2 ) + m(3) In 



E + ap + 13m 
E — ap + (5m 
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The limit \imF(x, A) = F(oo,A) can be calculated in an evident form 



F(oo,A) 



47r 2 \q\a 2 { 



\aA(VA 2 



+ m 2 + (3m) + m 2 [£(A, a 1: h) + f (A, a 2 , 62) 



/(A 2 + m 2 )(l-a 2 ) Q r- -\ 

+ I + /^WA 2 + m 2 j In 



VA 2 + m 2 -aA + (3m 
VA 2 + m 2 + aA + /3m 



Here 



£(A,a,6) 
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In 



aVA 2 + m 2 -m + AVa 2 - 1 



and 



/3 + «V« 2 + /^ 2 - 1 
01 = ; 5 , a 2 = 



am — VA 2 + m 2 

(3 - a^a 2 + (3 2 -l 



1 — a 2 

(a 2 + (3 2 ) (a(3 + V« 2 + /? 2 - 1) , 

"1 — ~ ~ J ^2 — 



1 - a 2 

(a 2 + /3 2 )(a/3 - V" 2 + P 2 ~ 1) 



2(1 -a 2 ) ' z 2(1 -a 2 ) 

The asymptotic of F(p, A) for large A can be written as 



^, . x A 2 A „. . ^ / m 2 . 

F p, A = c 2 — + ci- + c + / p + O — , 
m 2 m \ A 2 ' 



where 



f(p) = ~ 



Ait 2 \q\a 2 Jo 



00 g— 2 - E |p|p 



2pa - (£(1 - a 2 ) + m/3) In 



+ ap + /5m 
— ap + /5m 



c , Ci, c 2 are the constants: 



m 



c 2 = 



87r 2 |gla 2 



2a + (a 2 - l)ln 



1 + a 
1 — a 



(3m 2 



' Cl = 4^w{ 2a - ln T^}' 



(16) 



(17) 



c = c 2 + 



m 



Ck 



An 2 \^\a 2 \{l - a 2 ) 



(3 2 h In [ - ai - ya 2 - 1] 6 2 In [ - a 2 - - 1] 



where we used the notations (15), (16), and O (jp^ —> for A — > 00. Requirement of JF^ 
being finite for A — > 00 means that £ depends on A, and for large A its asymptotic has 
the form 

A 2 A' 



m 



£(A) = 2 c 2 — + Cl - +X + O — . 



m c 



m 



A 2 



Here the parameter \ is the renormalized value of £ . 

Thus, if we denote JF = (JF 03 , JF 13 , JF 23 ), then for A — > 00 we obtain the following result 

J?(x) = e Sign(x 3 ) g ( c + f(x 3 ) - X 
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The asymptotics of function f(p) (17) are of the form 

am 2 e -2m\ P \ 



M 
M 



p->°° 8|g|(7rm|p|) 3 / 2 (l + /3) 
c 2 Ci 



l+0(l/p)), 



p^o 2m 2 p 2 2m\p\ 



c + c 2 + O(p). 



Hence 
f(x) 



e Sign(x 3 ) q c 



X 



am 2 e -2 m \x 3 \ 



2 8|?|(7rm|x 3 |) 3 / 2 (l + /3) 



;i+o(i/x 3 )) , 



•F(x) = — eSign(x 3 )g 

£3^ 



c 2 



+ 



Cl 



2m 2 x 3 2m\x 3 \ 2 



+ ^ - c 2 + 0(x 3 ) 



We see, that on large distances the field generated by defect is of the same form as 
the field of the plane in classical ED [30, 31], and on the small distance x 3 it has non- 
classical behavior T ~ eqc 2 /2m 2 xl. Let us consider 3 cases: 1) q = (k, 0,0, r) = q^; 
2) q = (0, k, 0, t) = q^; 3) q = (k,k,0,t) = (/ 3 \ If g = g^, the defect generates pure 
electric field: H 1 = H 2 = H 3 = 0, E 1 = E 2 = 0, £ 3 = -iF^-ik). For g = g( 2 ) the field 
is pure magnetic one: E-y = E 2 = E 3 = 0, Hy = H 3 = 0, H 2 = J-y 3 . For g = g( 3 ) there 
are magnetic and electric fields: Ey = E 2 = Hy = H 3 = 0, E 3 = H 2 = JF 13 . 

To present the obtained results for all three cases in a compact form we introduce the 
following notations: Fi = E 3 for q = q^\ F 2 = H 2 for q = q (1 \ and F 3 = E 3 for q = q^ . 
Then near the plane the strengths of fields are 



Fi(x 3 ) 



fi 



x3->o m 2 x\ m\x 3 



+ 



9i 



, (-l) i+1 e X K 



where, 



/i = 
/ 2 = 



era 



8tt 2 «? 

em 2 
167T 2 o;2 

/3 = 



'1 + a^)Arctg(ai) — ai 



> 0i = ,_9_ 2 [Arctg(ai) - , 



1 — al)\n 



1 + a 2 



enm 



1 - a 2 

2 



37r 2 (4 + r 2 + A 2 )' 
Here we used the notations: 

4k 



- 2a, 



93 



, 92 



47T 2 «f 

ef3 2 m 2 
8tt 2 



ln „ — Q2 — 2^2 



Oi 2 



AenXm 2 



3vr 2 (4 + r 2 + A 2 ) 2 ' 
4/t „ Act, 



■ a 2 



4 - k 2 + r 2 + A 2 ' ~ z 4 + k 2 + r 2 + A 2 ' 
For the large distance asymptotic we have the following results: 



1,2. 



Fj(x 3 ) = Wj + 
£3^00 



l) t+ ex^ 



-2m|a;3 



(4m7r|a; 3 |) 3 / 2 



{l + 0(l/x 3 )} 
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The constants ■Uj^for i — 1,2 are written as 



Ui = 



era 

47T 2 <X 2 



1 - (-l)^,; 2 + h 



bf In [ - af - V(4 i} ) 2 -1 



fe=i 



4°) 2 - 1 



ec^m 

i + A' 



where 



(1) _ A + iai\jfii - «i 2 - 1 (i) _ (ft 2 - ai 2 )(«i/3i - «^//3i 2 - «i 2 - 1) 



1 + ai 2 



2(1 + «! 2 ) 



a 



A-iai^-ai 2 -! (1) (ft 2 - ai 2 )(ai/?i + - ai 2 - 1) 

' 2 " 2(1 + ax 2 ) 

(a 2 2 + & 2 )(a 2 ft + yV + ft 2 -!) 



1 + ai 2 



(2) ft + « 2 \/«2 2 +ft 2 -~ ( 2 ) 

) °1 



l-a 2 2 



(2) _ P2 ~ «2yW + & ~ 1 , (2) 
a 2 — ; o ) °2 



2(1 -a 2 2 ) 
(« 2 2 + (3 2 2 ){a 2 (3 2 - ^a 2 2 + (3 2 2 -l) 



l-a 2 2 

For u 3 , v 3 one obtains the following expressions 

^ 2 / 32A 2 



2(1 -a 2 2 ) 



enm 



u 3 



3- 



3tt 2 (4 + r 2 + A 2 ) V (4 + r 2 + A 2 ) 



, u 3 



4e/t 



r 2 + (2 + A) 2 



5 Fields generated by simplest defects 

The obtained results demonstrate the non-trivial dependence of the field E and H on 
the parameters k, A, r. In the main approximation E and H are linear functions of %, 
and E = H = 0, if x = 0, k = 0. Let us consider the simplest non-trivial case k ^ 0, 
X = A = r = 0. ' 

Asymptotics of E and H for large and small x 3 are the following. If q = q^\ the 
defect generates pure electric field E 3 



E% 



em 



x 3 ^o 8tt 2 u 2 

£/ 3 



1 + o; 2 )Arctg(o;) -uj 



m 2 x\ 



- 2 



(18) 



rr 3 -»oo 47T 2 C<J 



em 2 4/t 
[Arctg(cj) — o;J , u 



A - K 



2 ' 



For q = g( 2 ) the field is pure magnetic 

,2 



#9 



em 



o 167r 2 o/ 2 



'2m 

l + ^)ln— 



2a;' 



2™2 



(19) 
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em 2 

Ho 



x 3 ^oo 87T 2 U 2 ' 



l-u' 



. Ak 
> u = TTT2- 



A + K 2 



For q = q( 3 \ E\ = E 2 = H\ = H 3 = 0, and asymptotics of the fields E 3) H 2 are of the 
form 



E 3 PS H 2 ~ 

OO oo 



£3 

enm 2 

127T 2 ' 



In all cases the fields are constant at large distances from the plane. It is in agreement 
with well known results for homogeneous charge and current distributions on the plane 
in the classical electrodynamics [30, 31]. The magnitudes of the fields are defined by the 
factors m 2 /16n 2 , m 2 /8n 2 , m 2 /An 2 . They are of the order 10 17 cmr 2 -4 4 • 10 17 cmr 2 . It 
is more than the typical surface density 10 16 cmr 2 of atoms on the boundary of solids. 
The field E(x 3 ) is a monotonic function of x 3 . Therefore comparing the large- and small- 
distance asymptotics (18), (5), we can obtain an estimation for maximal distance d m from 
the plane for which (18) can be used: 

1 . . . , Arctg(c<j) — lu 

>2 + c{u), c(u) = 2- 



m 2 d 2 n uj — (1 + u; 2 )Arctg(a;) 

Since < c{uj) < 1, we have d m ~ 1,4 ■ 10 _10 cm 4- 1, 7 • 10~ 10 cm. The ranges of validity 

x 3 < d rn of asimptotics (19), (20) can be found analogously: 

d m ~ 2, • 10- I0 cm 4 2, 4 • IQ- I0 cm for (19), and d m ~ 1, 7 • 10- 10 cm for (20). 

Important feature of the fields generated by the considered defects for q = q^ , q^ 
is that they are singular by k — ±2. It means that these values of parameter k are the 
phase transition points, where in the case q = q^ the electrical field E 3 is changed in a 
sudden way, and in the case q = q^ the field H 2 become infinite. This phenomena seems 
to be not very surprising since it is similar to the known supercritical effects induced by 
perturbation of Dirac field by attractive 5 potential which causes diving of the ground 
state into the Dirac see by finite value of the coupling parameter [33, 34, 35]. 

The points k = ±2 are stable in respect to transformation k — > k' = A/k, for 
which uj{k) — > uj(k') = —u(k), E 3 (k) — > E 3 (k') = —E 3 (k), and uj'(k) — > uj'{k') = 
uj'(k), H 2 (k) — > H 2 (k') = H 2 (k). For q = q^ in each point of space the magnitude of 
field E 3 , considered as a functions of k is restricted: |-E(k)| < |-E(2)| < oo for all values 
of K. 

For the A = x — T — the short distance asymptotic is of the form E,H ~ 
const [l/(m 2 x 2 ) — 2]. In this case the relative correction of the next to leading term 
appears to be independent on parameter k describing specific properties of the defect on 
the plane. 
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6 Conclusions 



Suggested model describes an infinite plane with homogeneous charge and current dis- 
tributions in the framework of QED. Specific properties of the physical system are char- 
acterised by additional term Sdef (action of defect) combined with the usual action of 
QED into the full action of the model. Sdef was chosen on the basis of general princi- 
ples of QED: locality, gauge invariance and renormalizability of the theory. It contains 4 
parameters which can be fixed by normalisation conditions. 

The calculation of the leading order effects shows that mean field induced by Sdef has 
classical behaviour at large distance from the plane. Corresponding asymptotic can be 
used as normalisation conditions pinpointing interplay of parameters describing the fields 
of a plane in classical ED [30, 31] with ones of considered model. Along this way one can 
express four parameters of Sdef in terms of the effective charge and current densities and 
constants characterising macroscopic properties of material of the plane. 

On short distances the fields E, H are singular as functions of distance x from 
the defect : E,H ~ const/x 2 . Estimating energy density with usual classical formula 
W = (E 2 + H 2 )/2, one obtains its behaviour for (18), (19) and (20) at short distances 
as W(x) ~ c/x A + c'/x 2 , where c, c' are finite independent from cut-off constants. This 
singularities representing physical peculiarities of the model could be predicted with di- 
mensional analysis. It is similar with one found for scalar field under Dirichlet or Neumann 
boundary conditions on a single plate [36], akin local effects near surfaces can be observed 
in different geometries ( see, for example, [37, 38, 39]). Our model predicts dependence 
of c and d on parameters of the material plane. 

Properties of material film described by Sdef are defined by classical 4-current (given 
by ^(^4)) and particle density and 4-current of quantum of electron-positron field (given 
by S '\ q (ip , ip)) ■ Thus, it is assumed that the film consists of infinitely heavy charged 
particles, electrons and positrons uniformly filling the plane. Fluctuation of Dirac field in 
QED generates an effective width of the film. It is of the order 10~ 10 cm. In more realistic 
model one should describe the heavy charges in the framework of quantum fields theory. 
Simplest model of this kind could be QED with additionally proton-antiproton field. 
Fluctuations of this field generate anomalous electromagnetic field at a distance of order 
10~ 13 cm (because proton is 2000 times heavier then electron). Thus, they are insufficient 
for the effects at a distance of order 10~ 10 cm calculated in this paper. The leading 
singularity of the electromagnetic field near the plane is mass independent. Therefore in 
virtue of opposite charges of electron and proton, it is canceled for neutral defect with the 
same parameters A, q for both fermion fields, and at the distances x less than 10~ 13 cm 
their fluctuation generates the fields of the form < E >~ (M — m)C e /\x\ + (M 2 — m 2 )C' e , 
< H >~ (M — m)C h /\x\ + (M 2 — m 2 )C' h , where M (m) is the mass of proton (electron), 
C e , C' e , C h , C' h are functions of parameters A, q, and C e = C h = if A = 0. 

In the model of homogeneous defect the the surface energy density is infinite because 
of ultra-violet divergences, and needs a regularization with cut-off impulse A forming 
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effectively a discreet lattice structure with characteristic scale ~ A -1 in the coordinate 
space. It can be considered as a model of film consisting of atoms. Analysis of the function 
F(x 3 ,A) used for calculations of the fields J r ^(x 3 ) show that fields near the plane for 
discreet defect J r IJiU (x 3 ) can be described if we replace in results for continuous defect 
the singularities of ^^{x^ generated by fermion field of mass uif as follows: \jx\ — > 
f2(x 3 ,n) = (1 - (1 + nx)exp{-fi\x 3 \})/x%, l/\x 3 \ -> fi(x 3 ,n) = (1 - exp {-/j,\x 3 \})/\x 3 \ 
where \i = min{A,mj}. Such approximations are valid for distances smaller as ~ 
Thus, for neutral discreet defect with A^ 1 ~ 10~ 8 — 10~ 7 cm the behaviour of the fields 
^FpLu{x 3 ) on the distances ~ A -1 can be approximately described by main terms of short 
distance asymptotic of the homogeneous model. It is not surprising since there is no sharp 
cutoff dependence in the model of renormalizable quantum field theory. 

In our paper we have restricted ourself with the simple problem to calculate the mean 
electromagnetic field generated by perturbation of QED vacuum by infinite plane film. 
It is important to note that a finite physical observable is extracted for a system of just 
one isolated plane in distinction to ordinary Casimir effect. Obtained formulas predict 
dependence of Casimir-Polder forces on charge and current densities which can be proved 
experimentally. In the the critical point k = ±2 the magnetic field becomes infinite, and in 
the vicinity of critical point it must be observable near the defect surface with overcoming 
exponential suppression from masses of fermions. Other way to test the proposed model 
is the usage of obtained results for calculation of some physical feature of the film (for 
example, the spectrum of bound states of charged particle in anomalous electric potential 
of neutral defect) which can be investigated experimentally. In the framework of the 
model one can study scattering of electrons and photons on the defect and obtain results 
for experimental verifying. 

We hope that the proposed approach can provide a deeper insight into the nature of 
quantum phenomena of interaction of macroscopic bodies with QED fields. 
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